Experimental investigation of Wigner's reaction matrix for 
irregular graphs with absorption 

Oleh Hul^, Oleg Tymoshchuk^ , Szymon Bauch^, Peter M. Koch^, and Leszek Sirko^ 
^Institute of Physics, Polish Academy of Sciences, 
Aleja Lotnikow 32/46, 02-668 Warszawa, Poland 
"^Department of Physics and Astronomy, 
State University of New York, Stony Brook, NY 11794-3800 USA 
(Dated: September 22, 2005) 

Abstract 

We use tetrahedral microwave networks consisting of coaxial cables and attenuators connected by 
T-joints to make an experimental study of Wigner's reaction K matrix for irregular graphs in the 
presence of absorption. From measurements of the scattering matrix S for each realization of the 
microwave network we obtain distributions of the imaginary and real parts of K. Our experimental 
results are in good agreement with theoretical predictions. 
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Pauling introduced quantum graphs of connected one-dimensional wires almost seven 
decades ago l|. Kuhn used the same idea a decade later ^ to describe organic molecules by 
free electron models. Quantum graphs can be considered as idealizations of physical networks 
in the limit where the lengths of the wires greatly exceed their widths; this corresponds to 
assuming that propagating waves remain in a single transverse mode. Among the systems 
modeled W.ua„tu. g.aph. a.a a... electromagnetic optical waveguMe^ M 
systems [5|, 16|] , quantum wires P, |8| and excitation of fractons in fractal structures 
Recent work has shown that quan tum graphs provide an excellent system for studies of 



quantum chaos 
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211]. Quantum graphs with external 



15|. Quantum graphs with absorption, 



2lj, but 



leads (antennas) have been analyzed in detail in 

a more realistic but more complicated system, have been studied numerically in 
until now there have been no experimental studies of the effect of absorption. 

This paper presents results of our experimental study of distributions of Wigner's reaction 
matrix [22|] (often called in the literature just the K matrix 23|) for microwave networks 



that correspond to graphs with time reversal symmetry [jS = 1 sjTumetry class of random 



matrix theory 



24( 1) in the presence of absorption. For the case of an experiment having a 



single-channel antenna, the K matrix and scattering matrix S are related by 



(1) 



1 + iK 

The function Z = iK has direct physical meaning as the electrical impedance, which has 
been recently measured in a microwave cavity experiment |25| . For the one-channel case the 
S matrix can be parameterized as 

S = VRe'', (2) 

where R is the reflection coefficient and 6 is the phase. 

After seminal wotU of L6pe. Mello a.d Selig.a„ Q catne theotetical studies of the 

properties of statistical distributions of the S matrix with direct processes and imperfect 
coupling 28, 0, Is^. A recent experiment investigated the distribution of the 5* matrix 



for chaotic microwave cavities with absorption 31|- The distribution P{R) of the reflection 



coefficient R in Eq. (2), at the beginning investigated in the strong absorption limit |32| . 
has been recently known for any dimensionless absorption strength 7 = 27rr/A, where F 
is the absorption width and A is the mean level spacing. For systems with time reversal 



symmetry (/3 = 1) Mendez-Sanchez et al. |33| studied P{R) experimentally, and Savin et 



al. [34] found an exact formula for P{R)- For systems violating time reversal symmetry 
(/3 = 2), Beenakker and Brouwer Ssj] calculated P{R) for the case of a perfectly coupled, 
single-channel lead. 

In our experiment we simulate quantum graphs with microwave networks. The analogy 
between them is based on the Schrodinger equation for the former being equivalent to the 
telegraph equation for the latter We call them microwave graphs. Measurements of 



;he scattering matrix for them were stimulated by 



38l | and the pioneering measurements in 



37|. 



A simple microwave graph, the tetrahedral case, consists of six coaxial cables (bonds) that 
meet three- at-a-time at = 4 different T-joints (vertices). Each coaxial cable consists of an 
inner conductor with radius Vi separated from a concentric outer conductor with inner radius 
r2 by a homogeneous, non-magnetic material with dielectric constant e. The fundamental 
TEM mode that propagates (the so-called Lecher wave) down to zero frequency exists 
because the cross section of the cable is doubly connected [36|, p. 253]. For frequencies u 
jelow the onset of the TEn mode in a coaxial cable, which propagates above Uc — ^(^^^^^^)^ 
36(1, the cable is single mode: only the TEM mode propagates. For SMA-RG-402 coaxial 
cable, which has ri = 0.05 cm, r2 = 0.15 cm, and e ~ 2.08 (teflon dielectric), single-mode 
propagation occurs below 32.9 GHz. 

An (ideal) microwave graph with no aborption and no leads to the outside world is a closed 
(bound) system. The presence of absorption and/or leads to the outside world creates an 
open system. Because the coaxial cables are lossy, we may vary absorption in the microwave 
graphs by changing the length of cable(s), [21], by adding one or more (coaxial) microwave 
attenuators, or by changing the coupling to the outside world. 

Figure 1 shows our experimental setup for measurements of the single-channel scatter- 
ing matrix S for tetrahedral microwave graphs. We used a Hewlett-Packard model 8722D 
microwave vector network analyzer to measure the scattering matrix 5* of such graphs in 
two different frequency windows, viz., 3.5-7.5 GHz and 12-16 GHz. As the figure shows, at 
one of the four vertices we used a 4-joint connector to couple the microwave graph to the 
vector network analyzer via a single-channel lead realized with an HP model 85133-60017, 
low-loss, flexible microwave cable; the other three vertices consisted of T-joints. The plane 
of calibration in the measurements was at the entrance to the 4-joint connector. Note that 
the microwave graph in Fig. 1 has a microwave attenuator in one of its bonds. 
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FIG. 1: A diagram of experimental setup used to measure the scattering matrix S of tetrahedral 
microwave graphs with absorption. Absorption in the graphs was varied by changing the attenuator. 
The vector network analyzer used for all measurements was an HP model 8722D. 

To investigate the distributions of imaginary and real parts of the K matrix we measured 
the scattering matrix S for 184 different realizations of tetrahedral microwave graphs having 
a microwave (SMA) attenuator in one of the bonds. For each graph realization, which was 
obtained either by the replacement of the bonds or putting an attenuator to a different 
bond, the scattering matrix S was measured in 1601 equally spaced steps. The total optical 
lengths of the microwave graphs, including joints and the single attenuator, was 196.2 cm 
when a 3 dB, 6 dB, or 20 dB attenuator was used, whereas it was 197.4 cm when the 10 
dB attenuator was used. To avoid degeneracy of eigenvalues in the graphs, we chose optical 
lengths for the bonds that were not simply commensurable. 

Figure 2 shows the modulus \S\ and the phase 9 of the scattering matrix 5* of a tetrahedral 
graph with 7 = 3.6 (in panels (a) and (b), obtained with use of a 3 dB attenuator) and one 
with 7 = 6.8 (in panels (c) and (d), obtained with use of a 20 dB attenuator); both cases 
cover the frequency range 12-13.4 GHz. The lengths of corresponding bonds in the two 
graphs were the same. Direct processes are present in the scattering because the microwave 
vector analyzer was connected to the graphs by the 4-joint connector. For each, individual 
realization of the graph we may estimate them from the average value of the scattering matrix 
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FIG. 2: Panels (a) and (b) show, respectively, the modulus \S\ and the phase 9 of the scattering 
matrix S measured for the graph (see Fig. 1) with absorption parameter 7 = 3.6 (see the text) 
over the frequency range 12 - 13.4 GHz and with use of a 3 dB attenuator. Panels (c) and (d) 
show corresponding measurements for a graph with 7 = 6.8 over the same frequency range and 
with use of a 20 dB attenuator. The total optical length, 196.2 cm, of both microwave graphs was 
the same, including joints and the attenuator. 



{S). Our measurements averaged over all realizations of microwave graphs gave {S)av — 
0.47 ± 0.03 + i{-OM ± 0.04), where i = The experimental value for ^ 0.47 

is close to a theoretical estimate for the modulus of the vertex reflection amplitude 

IpI = 0.5 for a 4-joint connector with Neumann boundary conditions. 



P=--l, (3) 

where = 4 is the number of bonds meeting at the vertex in question. 

Equation (1) holds for systems with absorption but without direct processes. The case 



of imperfect^ coupling \{S)\ > and direct processes present can be mapped onto that of 
perfect one [23.] by making the following parametrization, 

9 - (A) 

- 1-\{S)\S' 

where Sq is the scattering matrix of a graph for the perfect-coupling case (no direct processes 
present). 

For systems with time reversal symmetry (/? = 1 in equations below), the distributions 
P{v) of the imaginary and P{u) of the real parts of the K matrix 23| are given by the 
following interpolation formulas: 

AT g— a 

P{v) = —^=-^{A[Ko{a) + Ki{a)]a + V^Be-^), (5) 



and 

where —v = ImK < and u = KeK are, respectively, the imaginary and real parts of 
the K matrix. The normalization constant is Np = a {AT{P/2 + l,a) + Be^°)~~^ , where 
a = 7/5/2, r(a;, a) = dtt^~^e~* is the upper incomplete Gamma function, A = — 1 and 
B = 1 + a — e". In Eq. (5) the variable a = ^{^/v + l/y/vY and Kq, Ki are MacDonald 



functions. In Eq. (6) D{z) = dq^l + z{q + ^-1)6-^^(9+9 and u = Vu^ + 1- 

Figure 3 shows experimental distributions P{v) for three mean values of the parameter 
7, viz., 3.8, 5.2, and 6.7. The distribution for 7 = 3.8 is obtained by averaging over 69 
realizations of microwave graphs having 7 within the window [3.5, 4.1]. The distribution for 
7 = 5.2 is obtained by averaging over 60 realizations of microwave graphs having 7 within the 
window [4.7, 5.6]. The distribution for 7 = 6.7 is obtained by averaging over 55 realizations 
of microwave graphs having 7 within the window [6.3, 7.1]. We estimated the experimental 
values of the 7 parameter by adjusting the theoretical mean reflection coefficient {R)th to 
the experimental one (Rq) = {SqSq), where 

{R)th= C dRRP{R). (7) 







We also applied the following interpolation formula 

G 

1- 

'R) 



3l| for the distribution P{R): 



PiR) = + ^(1 - Rf'\ (8) 
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FIG. 3: Experimental distribution P{v) of the imaginary part of the K matrix at different values 
of the mean absorption parameter: 7 = 3.8 (open circles), 7 = 5.2 (full circles) and 7 = 6.7 (open 
triangles), respectively. Each corresponding theoretical distribution P{v) evaluated from Eq. (5) 
is also shown: 7 = 3.8 (dotted line), 7 = 5.2 (solid line), and 7 = 6.7 (dashed line), respectively. 



We offer the following comment on the validity of the Eq. (8). We used it instead of 
exact formulas (12-14) recently presented in 3J], which may be used to find the distribution 
P{R), because Eq. (8) is sufficiently accurate (see Fig. 1 in [s^) while allowing for much 
faster numerical calculations. 

Figure 3 also presents for comparison with each experimental distribution P{v) (symbols) 
the corresponding numerical distribution (lines) evaluated from Eq. (5). We see that the 
experimental distribution P{v) at 7 = 3.8 and at 5.2 agree well with their theoretical 
counterparts. However, the comparison for 7 = 6.7 shows some discrepancies, particularly 
in the range 0.3 < v < 0.8. 

We may use measurements of the distribution P{u) of the real part of Wigner's reaction 
matrix for an imporant and natural consistency check on our determination of 7. Figure 4 
compares experimental and theoretical P{u) distributions at the aforementioned values of 7, 
viz., 3.8, 5.2, and 6.7. Though each case shows good overall agreement between experimental 
and theoretical results, for all three cases the middle (—0.25 < u < 0.25) of the theoretical 
distribution is slightly higher than its experimental counterpart. According to the definition 
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FIG. 4: Experimental distribution P{u) of the real part of the K matrix at different values of 
the mean absorption parameter: 7 = 3.8 (open circles), 7 = 5.2 (full circles) and 7 = 6.7 (open 
triangles), respectively. Each corresponding theoretical distribution P{u) evaluated from Eq. (6) 
is also shown: 7 = 3.8 (dotted line), 7 = 5.2 (solid line), and 7 = 6.7 (dashed line), respectively. 



of the K matrix (see Eq. (1)), such behavior of the experimental distribution P{u) suggests 
a deficit of small values of |ImS'o|. We do not yet know the origin of this deficit. 

Though there are the small discrepancies we have mentioned, the good overall agreement 
between experimental and theoretical results justifies a posteriori the procedure we have 
used to determine the experimental values of 7. 

The distributions P{v) and P{u) of imaginary and real parts of Wigner's reaction matrix 



may be also found using the alternative approach described in [25|, l26[ . In these papers the 



radiation impedance approach was developed and used to obtaining the distributions of real 
and imaginary parts of the normalized impedance 

Re Zc + i{lm - Im Z^) , . 

^ R^Z^ 

of a chaotic microwave cavity, where Zc{r) = Zo{l + Sc(r)) /(I — Sc(r)) is the cavity (radiation) 
impedance expressed by the cavity (radiation) scattering matrix Sc(r) and Zq is the charac- 
teristic impedance of the transmission line. The radiation impedance Zr is the impedance 
seen at the input of the coupling structure for the same coupling geometry, but with the 
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FIG. 5: Experimental distribution P{v) of the imaginary part of the K matrix at different values 
of the mean absorption parameter: 7 = 3.8 (open circles), 7 = 5.2 (full circles) and 7 = 6.7 



(open triangles), respectively, calculated using the radiation impedance approach [25j, |26t]. Each 
corresponding theoretical distribution P{v) evaluated from Eq. (5) is also shown: 7 = 3.8 (dotted 
line), 7 = 5.2 (solid line), and 7 = 6.7 (dashed line), respectively. 



sidewalk removed to infinity. This interesting approach is especially useful in the studies 
of microwave systems, in which, in general, both the system and radiation impedances are 
measurable. However, it is not obvious how to use in practice this approach in the case of 
quantum systems. 

We used this alternative approach to find distributions P{v) and P{u) of imaginary and 
real parts of Wigner's reaction matrix for irregular tetrahedral microwave graphs. Wigner's 
reaction matrix can be simply expressed by the normalized impedance K = —iZ. The 
radiation impedance was found experimentally by measuring in two different frequency 
windows, viz., 3.5-7.5 GHz and 12-16 GHz of the scattering matrix Sr of the 4-joint connector 
with three joints terminated by 50 Q terminators. 

In Fig. 5 and Fig. 6 we show the distributions P{v) and P{u) calculated using the radi- 
ation impedance approach {25, 26|. As in the case of the scattering matrix approach, the 
experimental distributions are obtained at three values of the parameter 7 = 3.8, 5.2 and 
6.7. Figure 5 shows that the distribution P{v) of the imaginary part of Wigner's reaction 
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FIG. 6: Experimental distribution P{u) of the real part of the K matrix at different values of the 
mean absorption parameter: 7 = 3.8 (open circles), 7 = 5.2 (full circles) and 7 = 6.7 (open trian- 



26l | . Each corresponding 



gles), respectively, calculated using the radiation impedance approach 
theoretical distribution P{u) evaluated from Eq. (6) is also shown: 7 = 3.8 (dotted line), 7 = 5.2 
(solid line), and 7 = 6.7 (dashed line), respectively. 



matrix for 7 = 5.2 is in good agreement with the theoretical prediction [23|]. However, for 
7 = 3.8 and 6.7 the theoretical results are slightly higher than the experimental ones, what is 
especially noticeable at the peaks of the distributions. The experimental distribution P{u) 
of the real part of Wigner's reaction matrix presented in Figure 6 at three values of the 
parameter 7 = 3.8, 5.2 and 6.7 displays a very good agreement with the theoretical result. 
The comparison of Figures 3 and 5 and Figures 4 and 6 show that the distributions P{v) 
and P{u) evaluated by means of the radiation impedance approach are at the peaks slightly 
higher than the ones obtained by the scattering matrix approach, what may suggest that 
the influence of the phase of 5*^ on the distributions is not negligible 261 ]. 

In summary, using the scattering matrix approach and the radiation impedance approach 
we have measured distributions P{v) and P{u) of imaginary and real parts of Wigner's reac- 
tion matrix for irregular tetrahedral microwave graphs consisting of SMA cables, connectors, 
and attenuators. Use of different attenuators allowed us to vary absorption in the graphs in 
a controlled, quantitative way. For the case of time reversal symmetry (/? = !), the exper- 
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imental results for P{v) and P{u) calculated for both approaches at the same three values 
of the mean parameter 7 are in good overall agreement with theoretical predictions. 
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